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A THEORETICAL STUDY OF THE EFFECT OF FORWARD SPEED ON THE FREE-SPACE 

SOUND-PRESSURE FIELD AROUND PROPELLERS ^ 


By I. li. Garrick and Charles E. Watkins 


SUMMARY 

The sound-pressure field of a rotating propeller in forward 
flight in free space is analyzed by replacing the normal-pressure 
distribution over the propeller associated with thrust and torque 
by a distribution of acoustic pressure doublets acting at the 
propeller disk and subject to uniform rectilinear motion. The 
basic element used to synthesize the field is the pressure field of 
a concentrated force moving uniformly at subsonic speeds, for 
which an expression generalizing one of Lamb’s for the fixed 
concentrated force is given. This result is presented both for 
the moving and for the fixed observer. The strength of the 
doublet distribution is related to the thrust and torque distri- 
bution and to its various Fourier coefficients in a convenient 
way. The sound field is expressed by integration over the 
propeller disk, and also by integration over an effective ring, and is 
given both for the near pressure field and, in a simpler form, for 
the far field. Known results for the zero-forward-speed case 
present themselves in the special case of Mach number M—0. 
Some illustrative examples are calculated and discussed. 

INTRODUCTION 

The rotating propeller is the source of an intense sound- 
pressure field which can be associated with the periodic 
reactions on the medium arising from the distifbution of 
pressure rotating along with the blades. This pressure 
distribution consists in part of a distribution due to thickness 
of the blades, whose resultant force in subsonic potential 
flow is zero, and in part of a distribution due to angle of 
attack and camber of the blades, whose integrated effect 
includes the induced drag and corresponds almost wholly to 
the thrust and torque distribution over the blade. Another 
source of propeller noise may be associated with flow separa- 
tion and with friction or shear due to the boundary layer; 
both effects lead to vorticity shed into the wake and hence 
the designation vortex noise. The vortex noise and the noise 
due to thickness (where wave drag is not a large factoi-) are, 
however, for actual propellers normally of a considerably 
smaller magnitude than the rotational sound due to torque 
and thrust; hence only the latter effect will be considered in 
the present work. 

A large number of investigators have studied various 
phases of the determination of the sound or noise field of 


rotating propellers. In addition to the references cited, a 
bibliography is included of representative work on this 
subject. A simplification that has frequently been made in 
propeller-noise investigations is to limit the considerations 
to the static or standing propeller. The work of Gutin in 
1936 (ref. 1) represents a development of this type that makes 
possible a satisfactory prediction of the amplitude of sound 
pressure due to thrust and torque of a propeller rotating on a 
stand in still air. Although Gutin’s theory is applicable to 
the near oscillating pressure field of the propeller, his results 
in reference 1 are limited to the determination of the funda- 
mental and the first few harmonics at a distance far from the 
propeller, that is, several diameters away. The determina^ 
tion of the near pressure field, however, has been of concern" 
both from structural and ph 3 'siological considerations. 
Hubbard and Regier (ref. 2) extended the application of 
Gutin’s work to describe the oscillating pressure field and to 
determine the amplitude of noise at points near the propeller, 
in some cases within a blade chord length from the tip. 
The\- investigated anahflicall^^ the effect of several of the 
parameters that enter in the theory and also gave compari- 
sons with experiment which were quite satisfactory. 

The existing, theoretical work has found useful application 
for static conditions and for conditions of low forward speed, 
for example, near take-off. For conditions of high forward 
speeds, however, many pertinent questions have arisen as to 
the possible effects of the forward speed on the oscillating 
pressure, field of the propeller. A few investigators have 
examined phases of this problem; Bryan, Hart, Shirokov, 
Blokhintzev, and particularly Kuessner and Billings (see 
bibliography) may be mentioned, but most of this work 
seems incomplete or difficult to appl,y. It appeared desirable, 
especially for applications, to j’e-examine the theoretical 
problem for the propeller at forward speed so as to have it 
arise as a straightforward generalization of existing work for 
the zero-forward-speed case. The purpose of the present 
report is therefore to extend the theory of Gutin (ref. 1) and 
the work of Hubbard and Regier (ref. 2) so as to include the 
effect of subsonic forward speed of the propeller on the near 
and the far oscillating sound-pressure field caused bj^ torque 
and thrust of the propeller. 


* Supersedes N ACA TN 3018, “A Theoretical Study of the Effect of Forward Speed on the Fiee-Space Sound-Pressure Field Around Propellers” by I. E. Garrick and C. E. Watkins, 1953. 
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This report includes the following material: (a) The sound- 
pressure field associated with a uniformlj^ moving concen- 
trated force is given, with details of the explicit development 
contained in the appendix. This result serves as a basic ele- 
ment for synthesizing the total pressure field . It is expressed 
both for the case of the observer or field point considered to be 
moving along rectilinearly and uniformly with the pi'opeller in 
free space and for the case of the observer considered fixed and 
the propeller in uniform flight, (b) The disturbance forces 
associated with the pressure distribution acting on the me- 
dium in the plane of the propeller disk are presented in the 
manner similar to that of Gutin. (c) The sound-pressure 
field resulting from the combined use of (a) and (b) is shown 
in the. form requiring integrations over the propeller disk 
and in the simpler form making use of an effective propeller 
radius. Approximations valid for the far field, which reduce 
to Gutin’s results for the case of zero forward speed, are also 
given, (d) Some numerical examples are calculated and 
discussed, (e) Remarks are made on the usefulness and 
limitations of tlie analysis. 
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t time 
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EXTENSION OF A FORMULA IN LAMB’S “HYDRODYNAMICS” FOR THE 
SOUND-PRESSURE FIELD OF A FIXED CONCENTRATED FORCE TO 
THAT OF A MOVING FORCE 


- -On-th e-basis of -acoustic-considerations of-the classical-hydro- 
d.ynamic equations. Lamb (ref. 3)^ gives the pressure at 
any field point x, y, z associated with an external periodic 
force Ae^‘ acting in the x-direction and concentrated' at the 
location xi, yi, Si as 


^ ^ e ^ 
4ir bx s 


(1) 


where 


,S = V(x-Xi)^-f (y — I/l)^+(2-2l)^ (2) 

The concentrated force maj- be regarded as stemming from 
the limit of a distribution of an increasing pressure differ- 
ence over a decreasing area, whose product in the limit of 
zero area is equal to the force. For an arbitrary concen- 
trated force of components P'x, Fy, F^ the result generalizes as 
Lamb indicates to 
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1 

1 d 0 i) 

, A 

Ldx s 





( 3 ). 


^ Dr. Nelson H. Ke.inp of the Research Department of the United Aircraft Corporation has 
in private correspondence kindly called the attention of the authors to an error of sign e.\isting 
in Lamb’s formula in his transition from equation (8) to equation (9) on page 502 of his “Hydro- 
dynamics.” This error was not corrected in NACA TN 3018 but is corrected herein, although 
it has no significance for any of the results of interest involving root-mean-squarc pressures. 
See, also, Kemp, Nelson H.: On the Pressure Field of a Uniformly Moving Concentrated 
Force. Jour. Acous. Soc. of America (Letters to the Editor), vol. 26, no. 3, May 1954, 
pp. 450-451. 
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This formula shows that the pressure field obtained from the 
concentrated force has the character of a doublet or a dipole 
in its periodic form it has played a central role in the develop- 
ment of the noise and pressure fields of propellers due to the 
propeller torque and thriist. However, the formula refers to 
a concentrated force and to coordinates fixed in space. To 
treat the case involving effects of forward speed, it is expedient 
to obtain the required extension of the preceding formula. 

The extension of equation (3) required for the case of a 
uniformly moving concentrated force is given in equation 
(4), which follows. Details for its derivation are supplied 
in the appendi.x. I^et the concentrated force have compo- 
nents Fi, Fy, F/, let it be moving uniformly with velocity V 
in the direction of the positive ai-axis; and let the coordinate 
system also move uniformly with the same velocity. 'J'hen, 
with X, y, z now denoting the field point referred to this co- 
ordinate system and with Xi, yi, Zi denoting the coordinates 
of the concentrated moving force 



V cj 


V cJ 




S 

' i>y 

s 

^d2 

S J 


( 4 ) 


where 

S=^Kx—Xly-\rP^[(y—yl)^+{z—z{)^ 

M(x—Xi) + S . . ^ 

r (5) 

Figure 1 illustrates the geometric meaning of the quantities 
s, S, and a. The force located at 0(xi,yi,zi) and the field 



' ■’ Since the acoustical pressure is p = —p the dipole represented' by equation (3) Is an 
acceleration doublet or the time rate of change of a flow doublet. 


point at Q(x,y,z) are both in uniform motion with velocity 
V in the a:-direction. The distance OQ=s. The influence 
at Q at time t stems from the action of the force when it 
was at position P, where P is obtained from the relation 
0PI0Q=M, the forward-speed Mach number. The distance 
PQ gives then directly the phase radius <r. The perpendicular 
dropped from 0 onto PQ determines QR, which is equal to 
the amplitude radius S'. ForM=0 both S and o- reduce to 
the ordinary radius s. Comparison of equations (3) and (4) 
shows that the effect of the forward speed leads to replace- 
ment of s by S' in tbe amplitude and by a in the phase. 

The force variation of main concern herein is that of a 
harmonically periodic force having components F^, Fy, F^ 
vaiying in time as Equation (4) becomes 



where ^=co/c. Equation (6) is the basic relation to be made 
use of in the subsequent anal 3 'sis. 

DI.STURBANCE FORCES IN THE PROPELLER PLANE 

Consider, as does Gutiii (ref. 1), a propeller disk Oriented 
so that the axis passing through the center of the disk coin- 
cides with the 2 i-axis and let the propeller be assumed to 
rotate, as in figure 2, in the 1 / 2 -plane (x=0) with positive 
values of x corresponding to points ahead of (and negative 
values of x corresponding to points behind) the propeller 
disk. T'he propeller is considered to move uniform!^" with 
velocity V in the positive a;-dircction. Points in the propeller 


2 



Figure 2. — Propeller disk and coordinate system. 
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disk are designated by ( 0 ,y,, 2 i) or, in polar coordinates 
(as shown in fig. 2), by 

yi=/'cos0| 

(7) 

2 i=r sin dj 

For definiteness, let the propeller be rotating counterclock- 
wise as seen by an observer looking into the propeller toward 
the slipstream. 

Each element of the propeller is acted on by the surface 
pressure distribution and this distribution may be resolved 
into a thrust force in the direction of the axis of rotation 
(the a;-axis) and into a force associated with the torque 
which acts about the axis of rotation opposing the rotation. 
Equal and opposite reaction forces to these are exerted on 
the medium. The points of application of these forces arc 
imagined to act in a single plane designated as the plane of 
rotation or as the propeller disk. 

The reaction of the surface pressure distribution of the 
rotating propeller on the medium at any instant is to be 
replaced by fixed periodic forces acting at the propeller disk. 
The propeller disk itself will be considered covered with the 
necessary singularities in the pressure or acoustic radiators 
of proper strength and harmonic content to correspond to 
this normal-pressure distribution. These singularities will 
be seen to be acceleration sources (acoustic radiator of zero 
-order)- for- those -forces- which - act- symmetrically- on -botli.-- 
sides of the blade — that is, whose net force over each element 
of the blade is zero — and to be acceleration doublets (acoustic 
radiator of order one) for those forces acting antis.ymmetri- 
call,y — that is, whose net force over each element of the blade 
corresponds to the difference in pressure over both sides of 
the blade. The doublet distribution is that which is needed 
to represent the thrust and torque distribution, in particular 
that part of the thrust and torque distribution associated 
with pressures acting normal to the blade surface. This 
part is practically^ all of it, except that arising from skin 
friction, (d'he effect of blade thickness may^ be taken into 
acepunt by introducing either flow sources or acceleration 
sources: the flow sources would lead to the sound-pressure 
field due to thickness by' tbe classical “piston” effect of the 
moving blade on the flow field; the acceleration sources 
would deal directly' with the assumed or known contribution 
to the pressure distribution over the blade due to thickness.) 

Consider an element of the propeller at distance r from 
the axis; let dr be its radial length and 6 its width measured 
in the projection onto the plane of rotation. Let the forces 
acting on the propeller element on each blade be A{r)dr in 
tbe a-xial flight direction and F{r)dr in the direction opposed 
to the direction of rotation. Equal reaction forces acting 
opposite to the flight direction and in the direction of rota- 
tion, resiiectively, are exerted on the medium. These quan- 
tities are related to thrust T and torque Q by- the relations 


(JT 

dr 


=BA (r) 




(8a) 


or 


T=B A{r) dr 

Q,=B r rF{r)dr 

Jo J 


(8b) 


where B is the number of blades and R is the length of each 
blade. 

The periodic impulse or reaction experienced at any* 
element of the disk may' be e-xpanded in a Fourier series. 
For simplicity' let the element considered first be located on 
the radial line 9=0 (the T/-axis) and afterwards be located 
arbitrarily'. To start with, assume that the forces are 
uniformly' distributed over the projected width 6 of the 
hlade element, that is, that the distribution of pressure 
difference over the blade chord is rectangular. Then in the 
area element r dr dd of the plane of rotation, the forces 

A(r)dr ^-^ and F{r) dr act on the medium during the 


time interval in which this element is eclipsed by the pro- 
jection of the propeller element. If the overlapping of the 

clement starts at t=0, it will end at t=T=--~ and the over- 

rQ 

lapping of the element will start again by' the next blade at 
2ir 

t=To=j^> where Q is tbe angular velocity' of the propeller. 

The rectangular- type forces experienced at the element of 
the disk-located-at-0=0-hy-its-periGdic- eclipse- by' -the-blade- 
may' be developed in a Fourier series: 


FAf)= 


U(r)'-fdr 

lo 


for 0<f<r 
for T<f<ro 


CD 


F^*it)=Ao+ Z] cos (mBQt—i,,;) 

m = l 


and similarly 


F-F(t)=^ 


( F{r) 


dd 


dr 


lo 


for 0<f<T 
for T<C_t<^ To 


( 9 ) 


where 


F2*(t)=Bo+'i2 B„, cos (mB^t—ri,,^ 

■m = l 

, 2 , , . r . rniTT , 

A „^= — /l(r) T sm dr dd 

rriT 0 To 

„ 2 r^, r . rriiTT , 

B,n=— (r) T sm dr dd ^ 

mw 0 To 

_ _mwT_mB^T 
Vrn o 

To ^ J 


(10) 


( 11 ) 


The constants and Bo which correspond to the instanta- 
neous average thrust and torque over the blade element and 
to the associated momentum shed into the slipstream do not 
give rise to sound and need not he expressed. The phase 
angles and which are small for the lower harmonics, are 
needed to preserve generality' of discussion since A„ cos 
corresponds to the coefficients of pure cosine terms and 
A,„ sin e„ corresponds to coefficients of pure sine terms 
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in the Fourier series. For the assumed rectangular distribu- 
tion, and r]^ may be observed to be equal to zero if the 
origin of time is chosen at the overlapping of the center line 
of the blade by the i/-axis since then only a pure cosine series 
suffices to express the distribution. 

The general formidation of interest where the disk element 
r dr do is located at an angle d may be expressed directly 
from these results, for, on a second area element r dr dO 
shifted with respect to the first b\' the angle 6 in the rotational 
direction, there act periodic forces of the same magnitude, 
but retarded by the time 0/U. The corresponding Fourier 
developments are 


F\{t)=Fi* (t—^\=Ao-\-'^A^ cos {mB9.t—mBd—t,n) 

\ IV m=\ 

F 2 {f)=Fj* (^—^=Bn-\S^,B„, cos imBQt—mBe—nm) 


r (12) 


The quantity 


mwT _rnBUt _mB b 
To 2 2r 


(13) 


is small for the lower harmonics, especially for blade ele- 
ments not near the hub. Blade elements near the hub, 
where r is small, are eliminated from consideration since 
they contribute very little to the air forces. One 

then replace in equation (11) sin FF!L1 so that, noting 

.To Tf» 

the relations in equations (8a), 

B 


To 


/I, 




A(:r)dr dd<=^~~ dr dd 
T 7T dr 

F(r)dr dr de 

7T Tr r dr 


(14) 


This appro.ximation, which weights the harmonic content 
of all the harmonics equally, may be observed to correspond 
to the case of the thrust and torque distributed over a zero 
blade width (that is, to the mathematical pulse sometimes 
termed the Dirac delta functioTi). This approximation 

FFUF. exceeds about 


T/'siny=0.707 
4 V 4 


; j=0.785^or' 


To 


becomes relatively less valid wheii 

when 4m— exceeds about unit 
of the 


To 


211 


harmonic exceeds 777 or. 
Bo 


that is, when the order 

roughly, 1/solidity. 

Moreover, the assumption of a uniform rectangular dis- 
tribution of the forces across the blade has been made for 
convenience, and other distributions may be treated if 
desired. vSome of the possible errors in the assumption 
leading to equation (14) are discussed briefly by Gutin in 
reference 1 and are shown to be generally negligible for the 
lower harmonics. Regier and Hubbard (ref. 4) also discuss 
this assumption in an illuminating manner. Figure 3, 
which is taken from reference 4, shows the relative harmonic 
content of different assumed distributions of the same total 
load for (a) the mathematical sharp pidse, (b) a triangular 



(a) Sharp pulse. 

(b) Triangular pulse. 77 ^- = 0.03. 

Z7rr 

(c) Rectangular pulse. ~^= 0 . 03 . 

(d) Rectangular pulse. ;t^ = 0.06. 

Figurk 3. — ^Effect of impulse shape on the relative amplitudes of the 

harmonics. 

‘ hat” pulse ;^=0. 03,. (c) a rectangular pulse =0.03, and 
Zirr lirr 

(d) a rectangular pulse 7^ =0.06. Blade widths of actual 

propellers tend to be between cases (c) and (d) for the most 
effective parts of the propeller disk. The sharp-pulse 
assumption generally tends to overestimate the magnitude 
of the higher harmonics. When the proper distribution is 
known, appropriate correction factors for the required 
harmonic may, of course, be applied to the magnitude of 
the results given by the pulse solution. 

The formulation in equation (12) need not be limited to 
the one in which the propeller force distribution is uniform 
throughout the propeller cycle. If interest should be 
attached to a nonuniform distribution, as occurs for the 
propeller yaived or pitched with respect to the flight path 
01' for one experiencing interference effects, it is readily 
possible to allow for these effects by permitting the ampli- 
tude of the distribution to become a function of 6. 

THE SOUND-PKE.SSURE FIELD 

The general expressions for the resolution of the forces 
associated with the thrust and torque on a radial blade 
element (for example, eqs. (12)) may be put in the usual 
convenient complex form (whose real part may correspond 
to the formulation of interest) ; 


W = 1 


(15) 


where w=mBS2 and where for the special case of the 
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rectangular-pulse distribution the coefficients are given as in 
equation (11). 

With the use of the approximations for A„ and given 
in equation (14), the periodic forces acting on the medium 
at the element r dr dd of the disk can be expressed for any 
given harmonic in the forms 




\_cIT 
TV dr 


dr de 


Fy sin e y 

7T r dr - 

F, „,=i i cos 0 dr dd 

7T /■ dr 


(16) 


In the derivation of equation (16), the phase angles e„, and 
rj„, which can generally be neglected for the lower har- 
monics, have been put equal to zero, a value which corres- 
ponds in the case considered (as has been mentioned) to the 
choice of 0=0 to correspond to the overlapping of the center 
line of one of the blades at the time t=Q. The index m indi- 
cates the order of the harmonic considered. As discussed in 
the preceding section, appropriate factors depending on the 
harmonic, number m, functions of the radius, and, for un- 
sjmimetrical loading, functions also of 6 ma,y (if known) be 
applied to the terms in equation (16). 

Pressure relations involving integrations over the propeller 
disk.— The pressure^at" any point iirthe fiee-s]iace field qiro- 
duced by these components of the uniformly moving peri- 
odic forces Fx, Fy, F^ is given in equation (6). With the 
use of the components of force given ih equation (16) and , 
by integration ovei- the propeller disk, the total oscillating 
pressure p for any given harmonic m (the index m will be 
hereinafter dropped)" is obtained from equation (6) as a sum 
of the pressure due to thrust pr and that due to torque p>Q 
as follows : ' 

P=Pt+Pq (17) 

where 


JL f" i 
IttOo Jo r 


'dT' 
dr * 


A iwt — rnBQ) 


d e' 
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1 (IQ 

dr ' 




dx S 
1 d e~’'‘ 


dr dd 


( 18 ) 


r bd S 


dr dd 


( 19 ) 


1 d 


and where, in the expression for Pq, the operation — — has 
been used for convenience to replace the equivalent operation 
The indicated differentiations with re- 


d d 

sm d ^ cos d -z^ 

oy. oz 

spect to X and d can be carried out by use of the expressions 
for O' and S given in equations (5) . 


d e~ 
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dr ® S 
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ikM ikx X 




drdd (21) 


and 


4Tr2jo Jo r dr 


• ’'“’/'ik B‘ 


^(ysind— 


z cos 9)dr (Id 

( 22 ) 


Another somewhat simpler expression for Pq can be obtained 
from equation (19) by integration by parts with respect to d 


pial !'R r 

Jo 


' 1 dQ 


0 P dr 


dr (Id 


S 


(23) 


Pressure relations involving effective ring approxima- 
tion. — Appreciable simplification may be achieved for calcula- 
tion purposes by making use of the approximations inherent 
in the assumption of an effective propeller radius Rg so that 
the integration with respect to r is avoided. Equations (21) 
and (23) then reduce to 


Vt= 


e’“' rjzfikM ikx :r\e 

“47?Jo \F^Ws'^~sy s 


and 


C'" Q ■ iM~' 






Vq I ^ ry 


S 


(Id (24) 


(25) 


where, in both S and a, the points yi and si have the values 
yi = Re cos d, Zi = Re sin d. In effect, equations (24) and (25) 
imply that the propeller disk has been replaced by an annular 
ring in which the entire thrust and torque are concentrated. 
■ d:'he effective'i’adius 'ok this ring “varies- somewhat, ’as' calcu- 
lation shows, with the load distribution and with the order 
of the harmonic. Doming (ref. 5) has shown by calculation 
in special cases for the static propeller that the ring approxi- 
mation is a reasonably good one. An effective radius of the 
ordei- of O.SR is considered reasonable for adoption in initial 
calculations. 

The magnitude of the root-mean-square pressure Pr„,s is 
of interest, since most sound-recording instruments arc 
calibrated in terms of it. The contribution of the harmonic 
of order m to the root-mean-squarc pressure is 

Vrms=’^\Pr + VQ\ 


or 


Pr„.s==^,{AMB'y>^- 

(mZf0-kfcv)+[T|(A/+|) 
— Q J sin {mBddrka) | 

sin i:mBd+ka)+^P^_(^M.+j 


where 


(26) 


-QSI cos QnBd+k<r)Y^ 


2_J I 

When rotational symmetry exists as for the condition of 
symmetrical loading, it is convenient in making the numerical 
calculations to let the field point be in the wy-plane, an 
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arrangement which can always be attained by suitable choice 
of the line 0=0. It may be desirable in some calculations to 
■divide the disk into several annular rings and to select 
effective radii for the various rings. 

Pressure relations for the far field. — A further simplifica- 
tion in the results can be reached if the distance s from the 
propeller disk {0,yi, z{) to the field point {x,y, 0) is large. Then 
from 

s‘^=x^+iy—yiy+zi^ 

there is obtained 

yyi ■ 


S«So — 


So 


(27) 


where So=^Jx^+y^ (the distance from the field point to the 
center of the propeller disk) . It also follows from 


that 


where 


S^=3^+'^^{y—yiy+P^Zi'^ 
o„ Q P^yyi 

O — 


(28) 




If terms of order IfS^ are neglected in comparison with 
terms of order 1 /S, the' souird pressure due to thrust is ob- 
tained from equation (24) by use of equations (5) and (28) : 




1 -.'.^(.so+Afx) 'iA: 




— imB6 + ik— ft- cos d 


From the known integral relation 


/: 


oiXCOS# 




do 

(29) 

(30) 


it follows that 


Vt- 


.-i^(So+ 


Mx) ik 


2wSo 

Similarly, from equation (25), 


P(m+|) /..(»•) (31) 






Observe that the argument of the Bessel function may be 
replaced by 

^^mBMrot f (33) 

oo ^0 

qr 

where Mrot= the Mach number corresponding to the 

c 

rotational speed at the effective radius. (The far-field 
approximation and the effective-radius approximation made 
use of in arriving at eqs. (31) and (32) need not, of course, 
be made simultaneously, since the required integration of the 
Bessel functions with respect to the radius may be carried 
out, at least numerically, without difficulty.) 

Introduce the fundamental frequency denoted by aii=/3n, 
so that 

u=mBQ=kc=7no:i 


The pressure magnitude for any harmonic m is then given by 


\p\ = \Pt+Pq\ 


mcoi 

27rCiS'o 



Be 

COii?/ 



(34) 


This result may be compared with that for the case of zero 
forward speed (M=0) given in reference 1 : 



. It can be noted that the forward-speed Mach number 
affects the torque term containing Q only in the replacing of 
So=V*^+P (the distance from the field point to the propeller 
hub) by the smaller distance So=-\lx^-\-^^y^. This s b- 
stitution occurs both in the argument of the Bessel function 
and in the outside factor 1/So- The thrust tei’m containing 
T is more strongl}'^ affected, since, in addition to this change, 
the field effect ahead of and behind the propeller disk is 

influenced by the term M-\—^ corresponding to a backward 

oo 

X ■ 

shift by -= =—M and also by an increase associated with the 

Oo 

factor l//3^. 


APPLICATION TO A SPECIFIC PROPELLER 

In order to give some indication of the effect of forward- 
speed Mach number on the sound pressure of a propeller, 
calculations based on equation (26) for the near field and on 
equation (34) for the far field have been made. For this 
purpose, a two-blade research propeller . having a 10-foot 
diameter and operating under the various conditions sum- 
marized in table I was chosen. As may be noted from this 
table, the propeller is assumed to operate at constant power; 
that is, the, power coefficient 6V and the torque Q ai’e held 
constant. As the forward-speed Mach number is changed, 
the blade angle and consequently the thrust T are changed 
according to the propeller charts of reference 6 so as to make 
the assumed conditions consistent with actual test operating 
conditions. Only the fundamental, the first harmonic m=l, 
is considered in these examples and the value chosen for the 
effective radius Eg is 0.8ft or 4 feet. It should be pointed 
out that the sound pressure computed from the data of table 
I would be obtained in pounds per square foot, since these 

TABLE I.— SEA-LEVEL OPERATING CONDITIONS FOR A 
10 -FOOT -DIAMETER TWO -BLADE PROPELLER 
(DATA FROM REFERENCE 6) 

[P=815 hp; Cp=0.10; 0=2,680 lb-ft,; fc = 0.29686] 


iM 

VhiD 

/Sr, deg 

Ct 

ij, percent 

T, ft-lb 

0 

0 

24 

0.11 

0 

1,850 

. 1 

.42 

26 

. 11 

46 

1,850 

.2 

.84 

28 

.094 

75 

1,600 

.3 

1.26 

33 

.065 

85 

1,090 

.4 

1.68 

37 

.052 

87 

875 

. 0 

2. 10 

42 

.041 

86 

690 

.6 

2. 52 

47 

.032 

80 

580 

. 7 

2. 94 

50 

.025 

70 

420 

.8 

3. 36 

53 

.0185 

62 

310 


3. 78 

56 

.0122 

46 

205 


327967—55 2 
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data are given in English units. These values have been 
converted to dynes per square centimeter in the results given 
in the figures multiplication b}^ the approximate conver- 
sion factor 480. 

Calculations for the near field. — Calculations of the root- 
mean-square pressure |p|^ based on equation (26) 

are made for various values of x in the range from — 0.5Z) to 
0.5Z1 along the line y=0.6D (6 ft), that is, along a line ^ of 
the propeller radius from the tip and extending a distance 
of 1 radius behind to 1 radius ahead of the plane of rotation. 
These results are shown plotted as a function of xjD for 
several Mach numbers in figure 4. It can be noted in this 
figure that, for each Mach number, two peaks generally 
appear, one ahead of and one behind the plane of rotation. 
For the set of conditions under consideration (see table I), 
the highest peak amplitudes of pressure occur just behind 
the propeller plane at values of xjD in the range from —0.15 
to —0.075. As the Mach number is increased from 0 to 0.4, 
the peak amplitudes decrease in magnitude but, as the Mach 
number is increased from 0.4 to 0.9, this trend is reversed 
and the peak amplitudes increase in magnitude, the peak for 
M=0.9 being about 1.4 times that for M=Q. Also, as M 
increases from 0 to 0.9, the point at which the highest peak 
pressure occurs moves somewhat nearer the propeller plane. 
Thus, the generally severe sound-pressure conditions at 
take-off (M « 0). ten d -to-be. alleviated .in flight_at . the lower- 
Mach numbers but ma}’' be reached again and even exceeded 
at the higher flight Mach numbers. The calculations for 
low Mach numbers presented in figure 4, it should be noted, 
are in substantial agreement with both calculated and 
measured results of reference 2. The high peak pressures 
obtained for the highest Mach numbers indicate that, with 
propeller-driven airplanes, the sound pressures generated 
near the tips of the propellers are of significance and of 
possible concern with regard to both structural considerations 
and passenger comfort. 

Calculations for the far field. — Calculations of the root- 
mean-square pressures based on equation (34) for various 
values of x along the line y=2D (20 ft) are shown plotted in 
figure 5. The trend with regard to Mach number in this 
figure is about the same as was noted in figure 4 for i/=0.6D, 
but the relative effect of the forward Mach number appears 
greater at a distance. However, because of the greater 
distances in figure 5 than in figure 4, the peak amplitudes 
are considerably less. 

Directional characteristics of the sound field. — For some 
purposes, especially with regard to calculations for the far 
field, it is desirable to consider the character of the sound 
pressure in terms of polar coordinates. For this purpose, 
substitutions can be made in equation (34) (and eq. (26)) 
as follows : 

x = So cos 4> y=So sin 4> (36) 


The quantities x/Sq and y/So in equation (34) (and eq. (26)) 
ma 5 '^ then be written as 


cos $ 


Sa -\/l— M^sin^$ 


sin $ 


(37) 


So 

With these substitutions, equation (34) becomes 




moil 


2Trcso-yJ\—M‘‘' sin^ 


|(mh 


cos 4> 


Vl— sin^ $ 


)- . 


Q 


Be j / kRe sin 4> \ 


(38) 


Calculations based on eqnation (38) for the Mach numbers 
of 0 and 0.8 for a constant value of So=2D are shown plotted 
as dashed curves in fignre 6. The solid curves in this figure 
represent results of calculations along the line y=2D 

obtained by replacing Sq by with y=2D in equation 

(30) for the various Alach numbers considered in figures 4 
and 5. A comparison of results of calculations along the 
line y—2D with those along the circle So=2D for M=0 and 
M—Q.S indicates that, for high Mach numbers (of the order 
of M=0.8), the peak pressures calculated along a line 
Constant -are -about -the same- as -those calculated along 
the circle with radius equal to the constant value of y. 
Observe the second pressure peak which has developed in 
the forward location at M=0.9. 

Separate components due to torque and to thrust of the 
sound field. — In order to give some indication as to the 
nature and proportion of sound pressure associated with 
each of the quantities, thrust T and torque Q, the root-mean- 
square pressures associated with each of these quantities for 
M=0 and M=0.8 are plotted for y=0.dD (6 ft) in figure 7. 
These plots show, as do equations (26) and (34), that the 
root-mean-square pressures associated with the propeller 
torque are symmetrically distributed with respect to the 
plane of the propeller for all Mach numbers, whereas the 
root-mean-square pressures associated with the thrust are 
symmetrical with respect to this plane only for M=0. For 
the particular propeller and operating conditions under con- 
sideration, the amplitudes of pressure associated with thrust 
for low Mach numbers are higher than those associated with 
torque, but for high Mach numbers the opposite is true. In 
the interpretation of figure 7, it should be recalled that the 
results depend on the assumed operating conditions and that 
the torque Q is the same (2,680 ft-lb) in both parts of the 
figure; whereas the thrust is 1,850 pounds at M=0 and 310 
pounds at M=0.8. Plots of this type can be used to obtain 
the sound pressure for various thrust and torque coefficients 
for a given propeller, since these coefficients appear as factors 
in equations (26) and (34) and hence can be normalized. 
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x/D 

Figure 4. — ^Distribution of the root-mean-square pressures for the fundamental (to= 1) of a two-blade, 10-foot-diameter propeller at several 
forward-speed Mach numbers. Cp = 0.10; fc = 0. 29686; t/=6 feet. (Operating conditions of propeller given in table I.) 
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Figure 
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440 


400 - 


M 

r,ib 

^ 0 

1,850 

.2 

1,600 

,4 

875 

.6 

580 

.8 

310 

.9 

205 


3601 


y-ZD 



. — Distribution of the root-mean-square pressures for the fundamental (»i=l) of a two-blade, 10-foot-diameter propeller at several 
fbrwmrd-speed Mach numbers. Cp = 0.10; i = 0.29686; y—20 feet. (Operating conditions of propeller given in table I.) 
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_ ^)_M=q^JT= 1850 1 _ (b) il/= 0.8: r=310 lb. 

Fk.iure 7. — Distribution of the root-incan-square total pressures and the root-mean-square pressures associated with thrust and torque for a 
two-blade, 10-foot-diameter propeller at forward-speed Mach numbers of 0 and 0.8. Cr> = 0.10; m=l; fc = 0. 29686; j/ = 6 feet. (Operating 
conditions of propeller given in table I.) 


It is of particular interest to note in the plot for M=0.8 
that the peak pressures associated with both thrust and 
torque are considerably greater than the peak pressure asso- 
ciated with the sum of these quantities. This result indi- 
cates that the phase relationship between the two components- 
is, in this case, such that each has a canceling effect on the 
other. Also in this case, the pressure associated with thrust 
may have its greatest value ahead of the plane of the pro- 
peller; whereas the total pressure has its greatest value behind 
the plane of the propeller. 

The discussioTi of the figures based on the specific examples 
illustrates the fact that the results do depend markedly on 
the assumed operating conditions. Moreover, only the 
fundamental harmonic {m=l) for the specific two-blade 
propeller has been illustrated. Trend studies on effects of 
higher harmonics, number of blades, and different operating 
conditions would be of considerable interest. Some pre- 
liminary calculations on the sound pressures associated with 
the fundamental of a six-blade propeller have shown a 
greater relative effect of the forward-speed Mach number. 

CONCLUDING REMARKS 

Expressions have been given for the sound-pressure field 
due to the distribution of thrust and torque for any given 
harmonic of a rotating propeller in uniform subsonic flight. 
The general expressions (eqs. (17) to (23)) involve integra- 
tions over the propeller disk; approximate expressions for 


the near fields (eqs. (24) to (26)) and for , the far field (cqs. 
(31),' (32), and 34)) involve integrations over a ring with an 
effective radius. The numerical examples have illustrated 
some free-space sound-pressure results for the fundamental 
of a specific two-blade propeller under various operating 
conditions at various forward speeds. 

It is pertinent to remark again on some of the limitations 
of the analysis, d'he analysis presented utilizes the torque 
and the thrust distributions, which may be given empirically 
or theoretically, in such a way as to require that they arise 
purely fi-om pressures acting normal to the surface of the 
blades. Empirical values of the thrust and torque include 
a conti’ibution, generally small, due to skin friction and to 
separated flow. Some caution is then needed in the use of 
the results. For example, the empirical torque term is 
somewhat lai’gcr than the toj-que due purely to normal 
pressures and, hence, the associated sound result due to 
torque may be overestimated slightly; similarly, the empiri- 
cal thrust term may slightly underestimate the sound due 
to thrust. At high tip speeds there may be significant con- 
tributions to the noise due to wave drag associated with the 
thickness. These contributions are mainly taken into ac- 
count b}'^ the effects of the wave drag on the torque. (Other 
sound effects of the thickness should be separately calculated 
and dncluded but these effects have not been explicitly 
presented herein.) 

Another way of looking at the approximation, and perhaps 
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a generally more convenient one for tlie study of trends, is to 
consider that the assumed torque and thrust distributions 
are actual theoretical oiies, obtained if necessary bj- adjust- 
ment of the blades of the propeller in a potential flow; and 
hence the sound-pressure field is that corresponding to the 
chosen thrust and torque in potential flow. Although the 
sound-pressure field and the aerodynamic velocity field have 
been considered as separate, it is of interest and of significance 
that the same concepts leading to the calculation of the sound 
field (made use of in the form of the acceleration potential) 
can lead to the linearized aerodynamics of the propeller in 
compressible flow, including the representation of the 
vorticity left behind in the wake as a result of spanwise 
variation of loading. The theoretical induced drag and the 
theoretical wave drag (which is a form of acoustic loss) are 
inherently included in the representation employed. 

It may be worthy of repetition that the pressure formulas 
as given tend to overestimate the contribution of the higher 
harmonics and that appropriate factors based on chordwise 
loading can be devised and ma\’ be required. Moreover, 
because of the theoretical (as well as the empirical) change in 


aerodynamic loading along the blade radially with increasing 
tip Mach number and forward-speed Mach number, the 
appropriate factors will change. In addition, the effective 
radius will be altered with the loading and with the harmonic 
considered. Thus, although 0.8 blade radius may be suitable 
as an effective radius for normal loading, a smaller value for 
the effective radius may be more suitable for conditions where 
unloading of the tip occui-s. The use of several effective 
radii corresponding to separation of disk into several rings 
may therefore be desirable. 

A calculation study of trends under different assumed 
conditions, the effect of unsymmetrical loading, or of dual 
propellers, effects of obstacles or boundaries on the free-space 
results, and experimental confirmations for the in-flight 
propeller are interesting matters for further investigation. 
It is also to be expected that effects of thickness will need 
to be taken into account in particular comparisons. 

Langley Aeron.^utical Laboratory, 

National .Advisory Committee for Aeronautics, 
Langley Field, Va., August 31, 1953. 



APPENDIX 


THE SOUND-PRESSURE FIELD ASSOCIATED WITH A UNIFORMLY MOVING FORCE 


DERIVATION. OF EQUATIONS (4) AND (6) 


Lainb (ref. 3, p. 502) derives from hydrod.ynamical princi- 
ples the differential equation satisfied by the pressure con- 
densation s* in the acoustic field associated with arbitrarily 
varying fixed forces acting on the medium 




\ da; dy bz ) 


(Al) 


where Lamb’s X, Y, Z are actually the extraneous forces per 
unit volume divided by the density p. (In effect, Lamb is 
dealing here with the acceleration potential.) In terms of 
the perturbation pressure j)=pc~s* the equation is 

□^p = divF 


dFx , bF„ . d/L 
dx ' by bz 


(A2) 


where F is the arbitrary force per unit volume having com- 
ponents Fx, Fy, Fx. Lamb^ shows that, if the jicriodic force 
is imagined concentrated on an infinitel.y small space 
at (xi, yi, 2 ,) and to be in the direction of x, the pressure at 
the field point (x, y, z) (the distance s from the location of 
the force) is given b.y 


where 


A A 

47t dx s 


(A3) 


. «=V(;«-a:i)^+ {y—y\Y+ {z-z\Y 


so that the concentrated force is equivalent to a double 
source or acceleration doublet whose axis is in the direction 
of the force. For a general harmonically periodic force F 
having components Fx, Fy, F^ varying in time as e'"' 



where k=~- 
c 


For an arbitrary time-dependent concentrated force F{t) 
located at (?,’), f), the pressure at field point (x, y, z) can be 
e.xqiresscd as 


p = 


-A cliv 

Itt 




47T 


dx 




F 

b 


by 


s ' bz s 


(AS) 


where the differentiations affect the result only through the 
variable s. 

The extension to a moving concentrated force may, within 
the framework of small-perturbation theory, be made in 
several ways — for example as mentioned in reference 7. One 
formal procedure utilized by Kiissner, for example, items 12 
and 13 of Bibliography, makes use of the invariance proper- 
ties of the wave equation and utilizes combined Galilean and 
Lorentz transformations. The following direct procedure 
given briefly in reference 7 is believed to be of intrinsic 
interest,- as it represents a simplification in Prandtl’s proce- 
dure (ref. 8) for the case of moving constant source distribu- 
tions, a procedure which consists of scheduling a succession 
of fixed soiirces in a path to act consecutively one after the 
other so as to represent in effect the desii’ed source moving 
along its path. Let the arbitrary concentrated force act 
only as an impulse during an infinitesimal interval at time 
t=T. The impulse ma.y be written as 

F(05(f-r) (A6) 


where the impulse function 6 ( t )=0 for tt^O and is charac- 
terized for f==0 as having unit area with respcct'to' r. The 
useful property of the impulse function of “sifting or select- 
ing” a value of a function is exhibited by the following rela- 
tion (see, for example, ref. 9, p. 61) : 



(A7) 


Let a succession of such impulses act, one following the 
other, in a path, points along which are given in space-fixed 
coordinates b.y . 

? = ?(r) J?=i7(r) f=f(r) 

The effect at time t of all such impulses which act before the 
time t is tlien given from Lamb’s result (eq. (A5)) as 





(A8) 


where, as has been defined, tbere is a nonzero contribution 
to the integral only for values of t defined bj' the character- 
istic relation 


i-r=l=Kl[x-^(T)YY[y-v{r)fY{z-t{r)f (A9) 

which expresses the distance between tbe source point and 
the field point in terms of the time of ti'avcl of the outgoing 
waves. 

The integral in equation (A8) corresponds to a summation 
of temporary fixed sources. To represent the case of uni- 
form j-ectilinear motion with velocity V in the positive 


2 Sue footnote 2, p. 2. 
14 
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2:-direction, let the sources be located on the J-axis and flow 
consecutively one after the other at the positions 


Vt r) = 0 


f=0 


(AlO) 


so that for t=— oo the source is located at J=— and for 
T = 0 the source is located at the origin. It follo\/s that 
the distance between source point and field point is 


,5 = V(3;-Fr)^-h7/ + 2^ . ■ 

Replace in. equation (A8) the variable r bj^ d where 


t-T — 


-e 


(All) 


(A12) 


By the sifting property of the 6-function (eq. (A7)) the 
integral will have a value only for 6=0 and equation (A8) 
becomes 


1 

“ dlV 
47T . 


P (r) dr-] 

[_ S dd^e=a 


(A 13) 


where the quantity within the bracket is now to be evaluated. 
Equations (All) and (A12) give a quadratic equation for 
T which, on choice of the solution that leads to and 
with 6=0, gives 


whore 


_ 1^ / Vx S 

V c 


8=^{x-vty+^^ (?/-+2^) 


(A14) 


/3=Vi-M' 


c 


Equations (All) and (A12) also give 
1 (It c 


s (Id es—V (x—V r) 

where replacing s bj- its characteristic value c{t—T) from 
equation (A9) and replacing t by its value given in equation 
(A14) yields, for 6=0, 


1 (It 1 

s (16~ S 


(A 15) 


where S is defined in equation (Al4). Equation (Al3) 
becomes 


= div 
47T 




(A16) 


This result is referred to space-fixed coordinates; it is 
finally desired also to express this result in terms of a field 
point {xo,yo,Z(^ of a coordinate system moving uniformly 
along with the source located at Xi, y-,, Zi so that one may put 

X—Vt = Xo—Xi y = yg — y, Z = Z„—Zi 

After this substitution is made, the zero subscripts may be 
dropped to yield as the end result corresponding to equation 


(A5) for the source and the coordinate system in uniform 
rectilinear motion : 


where 


I , '>H) 


■ s= V(x-x,)^+/3^ [(?/-?/.)"+ (2- Zif] 

M{x-Xi)-]-S 


(A17) 


iS=Vl-Af^ 

Equation (Al7) corresponds to the result given in equation 
(4) of the analysis. For the periodic harmonic force F 
varying as e‘“', it may be expressed in the form 


1 ,-r, e 


s 


(A18) 


which is equivalent to ecpiation (6) of the analysis. Equa- 
tions (A17) and (A18) are the sough t-for generalizations of 
equations (A5) and (A4) of the appendix, for the case where 
the disturbance and the field point are in uniform rectilinear 
motion. Comparison of these equations shows at once that 
s is replaced by S in the amplitude and by a in the phase. 
A geometric interpretation of these quantities is' shown in 
figure 1 and discussed in the analysis following equation (5). 

KEM.ARKS ON THE CASE OF THE MOVING DISTURBANCE AND THE 
FIXED OBSERVER 

A few remarks are in order on the significance of equation 
(A16). In this equation, the field point {x,y,z) is given in 
space-fixed coordinates while the disturbance force (or 
propeller) is in the in-flight condition moving with velocity 
V in the positive X-direction. Its location is given by 
Xi = Vt, yi = 0, 2] = 0 so that at t=0 it is located at the origin. 
It may be observed that the distance between disturbance 
and field point S is numerically the same whether given by 
the relation in equation (A16) or (A17). Hence, the pressure 
magnitude at any observer location for a sound-radiating 
element of the in-flight propeller is the same as that for the 
observer moving along with the propeller, provided the proper 
instantaneous distance between observer and propeller is 
used. However, there will be a difference in the frequency 
perceived by the observer. This frequency will be that for the 
case of the uniformly moving observer modified by the 
Doppler effect; thus, the frequency of each harmonic is modi- 
fied bv’ cl(c-\-Vr), where Hr is the component of the propeller 
forward speed (with proper sign: minus for approaching, 
plus for receding from the observer) in the acoustical direc- 
tion from the observer to, savq for the far approximation, 
the hub of the propeller. The acoustical direction is not 
quite that from the observer to the propeller location; it 
actually points from the observer to the location of the 
propeller when the sound which reaches the observer was 
emitted (direction QP rather than QO in fig. 1) . The Doppler 
frequency factor c/(c+H) is given geometrically in figure 1 

, . V QR S 

by the ratio 7^= — 
y X (T 
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